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FINAL  REPORT:  ARO  Grant  No.  DAAG29-77-G-0217, 
"Investigation  of  Load- Induced  Nonsynchronous 
Whirl  Instabilities  in  Rotating  Machinery" 


Background 

The  grant  for  this  investigation  was  made  to  the  Uhiversity  of  Florida 
(EIES)  in  the  amount  of  $70,469,  for  a period  of  two  years  beginning  in  Sep- 
tember 1977.  A third  year  of  funding  was  made  optional  at  the  discretion 
of  the  Army.  The  Principal  Investigators  were  Drs.  J.  M.  Vance  and  G.  N. 
Sandor. 

At  the  end  of  the  first  year,  Dr.  Vance  joined  the  Mechanical  Engineer- 
ing faculty  at  Texas  ASM  University,  leaving  the  grant  under  the  direction 
of  Dr.  Sandor.  Several  months  later,  the  grant  was  terminated  by  the  Army. 

This  report  gives  results  accomplished  up  to  the  termination  date, 
and  also  gives  certain  material  required  for  a possible  future  completion 
of  the  work. 

Objectives 

The  overall  objective  of  this  investigation  was  to  experimentally  verify 
(or  modify)  and  analytically  extend  the  two  available  theories  which  can 
explain  torque- load-dependent  nonsynchronous  shaft  whirling  in  turbomachinery, 
so  as  to  develop  an  improved  methodology  for  design  of  compressors  and  turbo- 
shaft engines. 

The  specific  technical  objectives  were: 

1.  To  experimentally  verify  in  a controlled  test  apparatus  the  math- 
ematical solution  of  reference  1,  called  "torquewhirl". 


2.  To  experimentally  verify  the  existence  of  the  aerodynamic  tip- 


clearance  effect  hypothesized  by  Alford  in  reference  7,  and  de- 


termine the  direction  of  the  induced  tangential  force. 


3.  To  extend  the  analysis  of  reference  1 to  include  the  effects  of  a 


more  general  type  of  shaft  flexibility. 


4.  To  derive  the  linearized  stiffness  and  damping  coefficients  re- 


sulting from  the  "torquewhirl"  generalized  forces,  for  use  in  a 


linearized  stability  analysis,  such  as  the  one  developed  by  Lund 


(reference  14)  for  multiple-disk  rotors  with  flexible  shafts  and 


other  realistic  complexities. 


5.  To  verify  the  usefulness  of  the  theory  by  incorporating  the  lin- 


earized "torquewhirl"  and  Alford  coefficients  into  a stability 


analysis  (such  as  the  computer  program  of  reference  14)  to  predict 


the  threshold  of  whirl  instability  for  a rotating  machine  with  an 


available  background  of  test  data  for  comparison  and/or  verification. 


Results  c-'d/or  Status  at  the  Termination  Date 


There  are  seven  tasks  described  in  the  grant's  Statement  of  Work  which 


are  designed  to  reach  the  above  objectives.  At  the  time  of  termination  of 


the  grant,  these  tasks  were  in  various  stages  of  completion,  as  described 


below.  Wherever  the  results  or  description  of  status  involve  extensive 


details,  they  are  relegated  to  referenced  appendices. 


Tasks  1 and  2:  Experimental  Verification  of  'Torquewhirl" 


A test  rig  capable  of  producing  torquewhirl  under  controlled  and  measured 


! 
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conditions  in  the  laboratory  has  been  designed.  The  basic  test  problem 
addressed  by  the  design  is  that  torquewhirl  normally  occurs  only  in  rotating 
machines  with  extremely  high  levels  of  load  torque  and  horsepower,  a condition 
which  is  expensive  to  simulate  in  the  laboratory.  This  problem  was  overcome 
by  designing  the  rig  with  very  low  danping  and  very  high  susceptibility  to 
torquewhirl. 

Drawings  are  now  ready  to  begin  shop  construction.  The  design  is  the 
result  of  an  extensive  optimization  study  to  produce  measurable  levels  of 
torque- excited  whirl  with  the  relatively  low  levels  of  torque  available  in 
the  laboratory  from  a 15  HP  aircraft  electric  motor.  The  design  includes 
a unique  concept  to  apply  the  torque  in  a controllable  manner.  See  Appendix 
A for  a detailed  description  of  the  design  procedure  and  the  resulting  test 
rig  design. 

Task  3:  Measurement  of  the  Alford  Force 

A test  rig  to  measure  Alfred's  force  has  been  redesigned  and  rebuilt 
with  a flexibly-mounted  fan  rotor,  adjustable  housing  offset,  variable- speed 
motor,  and  instrumentation  to  measure  the  deflection  produced  by  Alford's 
force.  The  instnmentation  includes  a special  filter  circuit  to  enable  the 
very  small  DC  deflection  signal  from  Alford's  force  to  be  measured  in  the 
presence  of  unavoidable  vibration  from  shaft  unbalance,  runout,  and  bearing 
roughness . 

Measurements  from  the  original  test  rig  did  produce  indications  that 
the  force  was  in  the  direction  predicted  by  Alford,  but  the  data  was  con- 
sidered inconclusive  because  of  surging,  large  standard  deviation,  and  the 


presence  of  a large  synchronous  vibration  component  in  the  signal. 

In  the  new  rig,  vibration  has  been  minimized  by  the  use  of  instrument 
quality  bearings,  and  by  precision  balancing  of  the  rotor. 

The  new  rig  was  just  ready  to  begin  providing  data  at  the  time  of 
writing  this  report. 

Details  of  the  rig  design  and  results  obtained  to  date  are  given  in 
Appendix  B. 

Task  4 : Comparison  of  Test  Data  with  Theory 

This  task  is  dependent  on  future  con?  let  ion  of  Tasks  2 and  3. 

Task  5:  Inclusion  of  a More  General  Shaft  Flexibility  in  the  "Torquewhirr’ 
theory 

This  task  is  a modification  of  the  original  Task  5,  in  keeping  with  a 
suggestion  by  one  of  the  Proposed  Reviewers  that  the  single- joint  shaft 
flexibility  model  in  reference  1 is  too  restrictive.  (Documented  by  letter 
correspondence  with  Dr.  Ed  Saibel). 

There  are  two  ways  one  can  take  a more  general  shaft  flexibility  model 
into  account.  Che  is  to  expand  the  three  degree  of  freedom  model  of  reference 
1 to  four  (or  five)  degrees  of  freedom,  allowing  the  shaft  to  bend  in  a 
general  manner,  and  attempting  to  find  an  exact  solution  to  satisfy  the  non- 
linear differential  equations  of  motion  (as  in  reference  1) . The  generalized 
coordinates  are  chosen  to  render  the  functional  solutions  tractable. 

The  second  approach  is  to  first  derive  generalized  forces  in  a cartesian 
coordinate  system  from  a virtual  work  expression  which  includes  the  four 
degrees  of  freedom  allowed  by  the  flexible  shaft.  Linearized  stiffness  and 
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damping  (Lund)  coefficients  are  then  derived  from  the  generalized  forces, 
and  can  be  used  in  a ccnputer  program  such  as  the  one  described  in  reference 
14. 

The  two  approaches  are  radically  different,  mathematically,  because 
of  the  different  choice  of  generalized  coordinates  to  describe  the  motion. 

The  former  approach  is  outlined  in  Appendix  C. 

The  latter  approach  is  the  subject  of  Task  6,  and  is  described  in  Ap- 
pendix D. 

Final  results  ham  not  been  obtained  as  of  the  termination  date,  but 
work  on  these  tasks  is  being  continued  at  Texas  ASM  University  by  Dr.  Vance 
and  Mr.  Brian  Murphy. 


Task  6 : Derivation  of  the  Linearized  Stiffness  Damping,  and  Cross -Coup ling 
for  Use  in  a Computer  Stability  Analysis  ~ 

These  coefficients  have  been  derived  for  the  'Torquewhirl"  model  of 

reference  1,  and  are  given  in  Appendix  "D". 

Derivation  of  the  coefficients  for  the  extended  flexible- shaft  model 

is  still  in  progress.  The  approach  is  outlined  in  Appendix  D. 


Task  7:  Load- Dependent  Stability  Predictions  for  a Machine  of  Realistic 
Complexity 

Part  of  this  task,  as  stated  in  the  Proposal . was  to  "improve  the  com- 
putational speed  and  accuracy  by  incorporating  the  method  of  reference  15 
to  generate  the  coefficients  of  the  characteristic  polynomial,  which  can 
then  be  numerically  solved  for  the  carplex  roots."  The  proposed  method 
involves  using  the  transfer  matrices  to  generate  the  polynominal  coefficients. 
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Since  this  is  conceptually  difficult  for  the  general  rotor  dynamic  problem 
with  four  boundary  conditions  and  complex  roots,  work  was  started  by  apply- 
ing the  method  to  torsional  vibration  problems.  Results  of  this  initial 
work  are  given  in  Appendix  D. 

The  major  part  of  Task  7,  a comparison  of  stability  predictions  with 
test  results  from  an  unstable  machine,  had  not  been  addressed  as  of  the 
terainatian  date,  since  it  is  dependent  on  results  from  Task  6. 
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Abstract 


Oh*  recently  published  Tbrquswhirl  theory  defines  a poaaihla 
aouraa  of  driving  anargy  which  nay  halp  cauaa  failura  dua  to  ncn- 
synchrcnous  shaft  whirling.  Xn  an  affdrt  to  investigate  this 
possible  driving  source,  a hast  rig  nuat  ba  daaignad.  The  daaign 
of  tha  taat  rig  ia  ftVMnwart  frtxa  tha  baginning  to  tha  pxaaant. 
Nuoaroue  problem  hava  baan  encountered  and  ara  presented  along 
with  thalr  propoead  solutions.  A datailad  dascrlption  of  tin 
test  rig,  as  now  anvisionad,  la  included. 
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conditions  in  the  laboratory  has  been  designed.  The  basic  test  problem 
addressed  by  the  design  is  that  torquewhirl  normally  occurs  only  in  rotating 
machines  with  extremely  high  levels  of  load  torque  and  horsepower,  a condition 
which  is  expensive  to  simulate  in  the  laboratory.  This  problem  was  overcome 
by  designing  the  rig  with  very  low  daiqping  aad  very  high  susceptibility  to 
torquewhirl. 

Drawings  are  now  ready  to  begin  shop  construction.  The  design  is  the 
result  of  an  extensive  optimization  study  to  produce  measurable  levels  of 
torque-excited  whirl  with  the  relatively  low  levels  of  torque  available  in 
the  laboratory  from  a 15  HP  aircraft  electric  motor.  The  design  includes 
a unique  concept  to  apply  the  torque  in  a controllable  manner.  See  Appendix 
A for  a detailed  description  of  the  design  procedure  and  the  resulting  test 
rig  design. 

Task  3:  Measurement  of  the  Alford  Force 

A test  rig  to  measure  Alfred's  force  has  been  redesigned  and  rebuilt 
with  a flexibly-mounted  fan  rotor,  adjustable  housing  offset,  variable -speed 
motor,  and  instrumentation  to  measure  the  deflection  produced  by  Alford's 
force.  The  instrumentation  includes  a special  filter  circuit  to  enable  the 
very  small  DC  deflection  signal  from  Alford's  force  to  be  measured  in  the 
presence  of  unavoidable  vibration  from  shaft  imbalance,  runout,  and  bearing 
roughness . 

Measurements  from  the  original  test  rig  did  produce  indications  that 
the  force  was  in  the  direction  predicted  by  Alford,  but  the  data  was  con- 
sidered inconclusive  because  of  surging,  large  standard  deviation,  and  the 
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presence  of  a large  synchronous  vibration  ccnponent  in  the  signal. 

In  the  new  rig,  vibration  has  been  minimized  by  the  use  of  instrument 
quality  bearings,  and  by  precision  balancing  of  the  rotor. 

The  new  rig  was  just  ready  to  begin  providing  data  at  the  time  of 
writing  this  report. 

Details  of  the  rig  design  and  results  obtained  to  date  are  given  in 
Appendix  B. 

Task  4:  Conparison  of  Test  Data  with  Theory 

This  task  is  dependent  on  future  conpletion  of  Tasks  2 and  3. 

Task  5:  Inclusion  of  a More  General  Shaft  Flexibility  in  the  'Torquewhirl" 
Theory 

This  task  is  a modification  of  the  original  Task  5,  in  keeping  with  a 
suggestion  by  one  of  the  Proposed  Reviewers  that  the  single- joint  shaft 
flexibility  model  in  reference  1 is  too  restrictive.  (Documented  by  letter 
correspondence  with  Dr.  Ed  Saibel) . 

There  are  two  ways  one  can  take  a more  general  shaft  flexibility  model 
into  account.  Che  is  to  expand  the  three  degree  of  freedom  model  of  reference 
1 to  four  (or  five)  degrees  of  freedom,  allowing  the  shaft  to  bend  in  a 
general  manner,  and  attempting  to  find  an  exact  solution  to  satisfy  the  non- 
linear differential  equations  of  motion  (as  in  reference  1) . The  generalized 
coordinates  are  chosen  to  render  the  functional  solutions  tractable. 

The  second  approach  is  to  first  derive  generalized  forces  in  a cartesian 
coordinate  system  from  a virtual  work  expression  which  includes  the  four 
degrees  of  freedom  allowed  by  the  flexible  shaft.  Linearized  stiffness  and 
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that  many  of  the  parameter*  used  In  the  design  of  a test  rig  are 
coupled,  lha  nature  of  the  coupling  is  rather  ocnplex  and  can  best 
be  seen  in  Figure  A-2.  There  are  further  couplings  between  par  one  tars 
that  are  not  readily  visible  from  Figure  A-2.  fbr  exanple,  it  la  not 
obvious  that  an  increase  In  the  moment  of  Inertia  about  the  point 
of  angular  deflection  results  in  an  increase  in  the  maunt  of  terse- 
pceier  required  to  produce  torquewhirl.  Nor  is  it  obvious  that 
this  seme  increase  in  the  moment  of  inertia  reduces  the  shaft 
esquired  for  torquewhirl.  Hidden  relations  such  as  these  make  it 
extremaly  difficult  to  get  a rig  design. 


U 
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I - aoaent  of  inertia  about  axis  of  rotation 
I • ament  of  inertia  about  point  of  angular  deflection 
Kq  - stiffness  refered  to  point  of  angular  deflection 
f - whirling  speed  ratio  \ 

Ws  ■ shaft  speed 
1 • shaft  length 

H P-  horse  power  required  to  produce  Torquewhirl 
if  • dating  ratio 


Design  Concepts 

Before  the  ocmplax  nature  of  the  par  one  tar  relations  was 
understood , a test  rig  wee  constructed  fay  another  investigative 
team.  Not  surprisingly,  this  rig  failed  to  produce  torquawhirl.  The 
reason,  of  aourse,  was  that  many  of  the  laportant  relations  were 
completely  overlooked.  A sketch  of  this  system  Is  shown  in  Figure  A-3. 
Vtiile  this  rig  was  basically  a failure,  it  did  provide  this  effort 
with  a starting  plaoa. 


It  was  known  beforehand  that  the  maximun  power  available  would 
be  around  15  H.P.  This  fact  put  limitations  on  the  system  design 
that  were  not  fully  understood  at  the  time.  The  first  attsnpts  at 
the  system  development  were  baaed  on  power  dissipation.  The  basic 
design  was  to  have  acme  sort  of  power  dissipation  unit  suspended  on 
a shaft.  This  shaft  was  to  be  coupled  to  the  main  drive  shaft  at 
the  point  of  angular  deflection.  (See  Fig.A-4)  The  outside  of  the 
dissipation  unit  was  not  to  turn  but  was  to  be  held  fay  springs  that 
were  grounded  to  the  support  frene.  The  reason  for  this  was  so  that 
the  whirl  motion  oould  be  ovserved  without  any  rotation.  This 
also  was  the  easiest  way  of  applying  a reaction  load  torque  to  the 
dissipation  unit.  The  type  of  dissipation  oould  be  either  aoulcmb 
friction  or  viscous  friction.  Vlsoous  friction  was  chosen  for  this  system. 
It  was  thought  that  an  Inpellar  oould  be  placed  inside  of  a housing 
that  was  supported  fay  the  sans  shaft  that  was  used  to  drive  the 
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impeller.  The  power  was  to  be  dissipated  by  puiping  the  fluid  from 
the  center  out  to  the  edges  and  then  recirculating  it  back  to  the 
canter.  This  method  was  dropped  because  of  the  many  problems 
associated  with  power  dissipation  calculations  and  the  possibility 
of  pup  stall  or  cavitation.  Other  fluid  puiping  ideas  were  abandoned 
for  the  sane  general  reasons.  The  second  approach  was  based  on 
viscous  shear.  In  this  idea,  the  housing  fit  very  close  to  the 
rotating  inpeller.  The  power  wes  to  be  dissipated  by  the  shearing 
of  a thin  film  of  fluid  between  the  housing  and  the  impeller. 

This  method  would  dissipate  s good  deal  of  horsepower,  but  required 
a large  diameter  impeller  because  it  was  desired  to  dissipate  approx- 
imately 15  H.P.  An  operational  speed  of  250  rpm  wes  arbitrarily 
assigned.  These  limitations  were  placed  an  the  systen  before  the 
relations  between  then  wes  known.  Because  of  the  low  rpm  of  the  system, 
a large  inpeller  wes  required.  Fbr  this  reason  the  idea  was  rejected. 
This  idea  may  be  more  viable  when  viewed  in  the  light  of  higher  speeds, 
which  appear  to  be  possible.  A review  of  this  idea  is  warranted  before 
the  go  ahead  is  given  to  any  other  design. 

Another  idea  based  on  viscous  friction  was  put  forward  at  a 
design  meeting  on  January  17,  1978.  In  this  idea  the  inpeller  took 
the  form  of  a conical  plug.  (See  fig.  A-5)  The  plug  fit  inside  of  a 
mating  housing.  The  shape  of  the  outside  of  the  housing  was  not 
determined  st  the  time  but  was  assumed  to  be  boxlike.  Power  dissipation 
in  this  system  wes  adjustable.  The  adjustment  wes  made  by  moving 
the  plug  in  or  out  slightly  and  thereby  changing  the  thickness  of 
the  visaous  shear  film.  This  is  basically  the  form  of  the  system  today. 
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At  the  timer  the  size  of  the  system  was  still  unknown.  In  an 
effort  to  get  sane  feel  for  hew  large  car  small  the  system  should 
be,  several  mathematical  models  were  set  ip.  An  equation  was 
developed  which  related  horsepower  dissipation  to  size,  film 
thickness,  viscosity,  shaft  speed  and  angle  of  plug  cone. 

Biis  equation  was  used  to  develop  horsepower  dissipation  vs. 
speed  curves  for  the  various  models.  1b  do  this,  the  visaosity  was 
assumed  to  be  56  micro  reyns  (Dow  Coming  silicon  fluid  1000) . The 
film  thickness  was  taken  as  .003  inches.  The  cone  half  angle  was 
taken  as  30°.  It  is  important  to  note  that  during  the  design 
process  the  value  of  f , the  Whirling  frequency  ratio  was  limited  to 
be  between  .4  and  .7  so  as  to  resonble  actual  field  date.  The  curves 
were  reviewed  in  light  of  this  limitation  on  f and  the  assured  horse- 
power available.  From  the  review,  one  system  was  chosen  as  moat 
desirable.  At  the  time  this  selection  was  made,  most  of  the  relations 
between  parameters  were  still  to  be  uncovered.  (Ftar  example,  though 
the  system  selected  from  the  curves  would  dissipate  an  adequate  amount 
of  horsepower  at  a certain  rpm,  it  would  not  be  dissipating  enough 
at  that  speed  for  torquewhirl  to  be  produced.)  The  selected  system 
was  then  analyzed  more  carefully.  The  effect  of  changing  the  shaft 
length,  the  effect  of  changing  the  speed,  and  most  other  parameters 
as  well,  were  investigated.  This  is  when  the  relationships  began  to 
appear  between  power,  speed  and  moment  of  inertia  about  the  coupling 
(located  at  a point  of  angular  deflection) . The  selected  system  was 
continually  modified.  The  purpose  of  all  modifications  was  to  increase 


* Note:  than  are  two  powers  with  which,  we  are  oonoamed.  One 
is  tha  power  tha  unit  will  dissipate  at  a given  rpn,  and  the  other 
is  the  power  required  for  torquswhirl. 


Problems  with  the  Basic  Unit 


Par  any  machine,  the  percent  of  critical  damping  that  will  be 
present  in  the  final  construction  is  very  difficult  to  determine. 

Ft>r  the  torquewhirl  test  rig  the  percent  of  damping  is  a very 
important  parameter.  It  has  a direct  effect  on  the  amount  of  horse- 
power that  is  required  for  torquewhirl.  Before  any  finalized 
design  can  be  developed  the  value  of  ^ must  be  determined.  To  do 
this  a rig  model  was  constructed  and  the  percent  critical  damping 
was  determined  by  the  log  decrement  method.  The  value  obtained  wes 
^ » .241.  This  very  low  value  enables  the  system  to 
function  at  much  higher  speeds  since  it  requires  less  power  to  produce 
torquewhirl. 


The  system  in  its  present  form  will  dissipate  approximately  12  H.P. 
This  brings  vp  an  obvious  problem.  Hew  to  get  rid  of  12  H.P.  worth 
of  heat.  It  was  at  first  hoped  that  the  thermal  capacity  of  the 
alunin’m  plug  and  housing  would  be  sufficient  to  absorb  the  heat  while 
the  necessary  readings  were  taken.  This  did  not  prove  to  be  the  case. 
Based  on  thermal  capacitance  alone,  the  operating  time  before  a 
temperature  in  excess  of  250°  F was  reached,  was  only  59  seconds. 

Seme  method  of  heat  transfer  would  have  to  be  devised.  Calculation 
for  free  convection  showed  no  improvement,  so  forced  convection  wes 
considered.  It  was  decided  that  fins  would  be  necessary  to  get  the 
rate  of  heat  transfer  required.  The  fins  are  to  be  machined  into  the 
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previously  cylindrical  housing  as  shown  in  Figure  A-7.  Forced  flow 
was  assuned  at  six  feet  per  second,  based  on  experimental  results 
from  a standard  blower.  Analysis  for  the  fins  was  carried  out  apart 
frcm  the  rest  of  the  housing.  Baaed  on  the  flow  velocity  and  an 
assured  steady  state  tanpazature  difference  of  175°,  a aonvective 
heat  transfer  coefficient  of  197.5  BfKJ/hr-ft^-°F  was  obtained.  The 
non-finned  portion  of  the  housing  was  > analyzed  based  on  the  flat 
plate  assumptions.  It  turned  out  that  eight  fins,  3/8  of  an  inch 
thick  would  be  sufficient  when  combined  with  other  surface  transfer 
rates.  The  forced  convection  is  to  be  developed  by  attaching  a snail 
fan  to  the  portion  of  the  shaft  that  is  rigidly  fixed  in  the  bearings 
above  the  coupling.  The  fan  will  turn  at  the  shaft  speed  of  the 
syston.  Air  will  be  drawn  ip  through  a large  tube  that  surrounds 
the  dissipation  unit  (aea  Fig.  A-8) . The  flow  velocity  can  be  controlled 
by  ports  at  the  bottom  of  the  tube.  Partitions  will  be  placed  in  the 
tube  to  ensure  that  the  flow  is  turbulent  before  it  reaches  the 
dissipation  unit. 

Another  problem,  not  quite  as  obvious  as  the  heating  problem,  is 
a result  of  the  centrifugal  effect  the  rotating  plug  Imparts  to  the 
viscous  film  layer.  This  centrifugal  effect  causes  the  fluid  to, 
in  effect,  be  pimped  out  from  under  the  plug.  The  net  result  is 
that  the  plug  is  sucked  down  with  more  and  more  force.  This  continues 
until  the  plug  seizes.  TO  alleviate  this  problem,  the  plug  is  fitted 
with  two  return  ducts  which  oan  be  seen  in  Figure  a-9.  These  ducts 
insure  that  the  fluid  will  be  able,  to  get  back  down  under  the  plug. 

The  angled  portions  of  the  ducts  act  as  centrifugal  putpa  In  ttamselves. 
The  use  of  this  solution  to  the  seizure  problem  requires  a level 
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of  fluid  to  cover  the  top  of  the  plug.  In  order  to  prevent  this  fluid 
_ front  being  slung  out  it  is  necessary  to  put  a cap  over  the  entire 

housing,  This  plate  serves  two  purposes.  First,  of  course,  it  keeps 
~ the  fluid  front  being  slung  out,  and  seoond,  it  provides  a seat  far 

the  bop  bearing.  This  top  bearing  is  included  to  ensure  alignment 
and  prevent  housing  wobble.  (See  Fig.  A- 10) 

In  order  to  make  adjustments  to  the  power  dissipation  it  is 
necessary  to  be  able  to  adjust  the  film  thickness.  In  the  original 
basic  design  this  vats  done  by  tightening  s nut  that  turned  with 
the  bottom  of  the  shaft.  Naturally  a more  desirable  method  would 
_ be  one  that  allowed  adjustments  while  the  rig  vats  in  operation,  fbr 

this  purpose,  a device  wes  developed  that  utilizes  two  wedges.  (See 
Fig.  A- 10) . The  top  wedge  is  fixed  with  respect  to  the  cover  plate. 

The  bottom  wedge  is  free  to  slide.  The  bottom  wedge  applies  pressure 
to  the  bearing  in  the  cover  plate.  This  bearing  sits  on  a shoulder 
- on  the  shaft.  It  has  a sliding  fit  between  its  outer  race  and  the 

cover  plate.  Nhan  the  lower  wedge  is  slid  in,  the  bearing  and  therefor 
the  shaft  and  plug  are  forced  down.  The  film  clearance  is  reduced. 

The  lower  wedge  is  controlled  by  tNO  thumb  screws,  one  on  either  aide, 
that  are  grounded  to  the  plate.  The  vedgea  have  an  angle  of  1.78°. 

— It  requires  1/2  inch  of  linear  wedge  motion  to  produce  .015  in.  of 

axial  plug  motion.  This  corresponds  to  full  adjustment  of  the  film 

r 

thickness.  TWo  bevel  springs  are  located  beneath  the  bearing  and 
_ contact  the  bearings  outer  race  and  the  plate.  The  purpose  of  these 

two  springe  is  to  aid  gravity  in  holding  the  film  gap  at  the  desired 

L 
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A good  deal  of  Marching  mu  required  to  find  a coupling  that 
Mould  handle  the  anticipated  68  ft.  lb.  of  torque  and  alao  the 
expected  4°  of  angular  miaalignnent.  Parallel  misalignment  features 

a 

found  on  most  aouplings  are  not  desirable  for  this  application.  The  * 

first  choioe  was  a universal  joint.  .However,  universale  do  not  transmit 
a constant  velocity  and  were  therefor  deemed  unsatisfactory.  (Amorous 
other  ootpling  styles  ware  investigated.  The  final  choioe  was  a 
Thomas  Modal  162.  The  Thomas  coupling  is  of  the  standard  disc  pack 
construction  and  will  handle  the  expected  loads.  One  possible  problem 

\ 

I 

that  requires  testing  is  the  added  danping  due  to  this  type  of 
ocupling.  It  is  hoped  it  is  not  sufficient  to  raise  the  horsepower 
requirement  out  of  the  range  of  the  selected  drive  motor.  Should 
this  danping  prove  too  large,  the  coupling  search  will  continue. 


for  powering  the  system,  a 14  H.P.  aircraft  generator  was  selected 
The  basis  for  this  selection  was  primarily  price,  but  also  because 
some  experience  has  been  had  with  this  type  of  power  drive  from 
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Summary 

The  primary  design  of  tha  tnrquawhirl  tost  rig  la  vary  naar 
oanpletion.  The  design  calculations  are  in  the  process  of  being 


checked  out.  ttaterlal  selection  has,  for  the  most  part,  bean 
completed.  Almost  certainly  there  will  be  inexpected  problem  that 
will  arise  during  and  after  construction.  These  will  be  dealt  with 
as  they  develop.  Every  precaution  has  bean  taken  so  that  nothing 
has  bean  overlooked. 


At  this  time  the  project  is  nearing  the  construction  phase.  What 
remains  to  be  dona  is  a final  check  of  all  calculations.  A recheck 
of  the  speed  and  power  requirements  for  the  system  with  all 
the  changes  included  in  the  calculations  and  the  neoasaary  parts  and 
assembly  drawings  will  be  done.  If  all  goes  as  expected,  construction 
will  begin  within  the  next  two  months.  Figure  A-10  is  a sketch  of  the 
system  as  now  designed. 
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Final  Design  of  Torquewhirl  Test  Rig 

The  design  of  the  Torquewhirl  test  rig  proved  to  be  far  more  complicated 
than  originally  assumed.  A thorough  examination  of  the  Torquewhirl  equations 
reveals  a high  degree  of  cross-coupling  between  the  variable  design  parameters. 
This  coupling  became  a major  problem  in  the  early  stages  of  the  design,  when 
initial  sizing  of  the  rig  was  necessary.  The  diagram  on  the  following  page 
(Fig.  A-1A)  gives  an  overview  of  the  linking  of  variable  design  parameters. 

Once  the  nature  of  the  parameter  coupling  was  understood,  the  design  process 
advanced  more  rapidly.  Thus,  changes  and  new  ideas  could  be  considered  in 
light  of  their  overall  system  effect.  (See  Appendix  A,  page  A3.) 

It  was  known  beforehand  that  the  maximum  power  available  would  be  around 
IS  h.p.  This  fact  put  limitations  on  the  system  design.  The  exact  nature  of 
these  limitations  was  not  fully  understood  at  the  early  stage  of  design.  The 
first  attempts  at  system  development  were  based  on  power  dissipation  alone. 

The  basic  system  was  to  have  some  sort  of  power  dissipation  unit  suspended  on 
a shaft.  This  shaft  was  to  be  coupled  to  a main  drive  shaft  at  a point  of 
concentrated  flexural  elasticity  (see  Fig.  A-2A).  The  outside  of  the  dissipation 
unit  was  not  to  rotate,  but  was  to  be  held  by  springs  which  were  grounded  to 
the  support  frame.  These  springs  afford  a method  for  measuring  the  torque 
developed  by  the  power  dissipation  device.  This  is  accomplished  by  observing 
the  rotational  displacement  of  the  housing  of  the  dissipation  device.  With 
this  type  of  setup  it  would  be  possible  to  observe  the  expected  whirl  motion 
without  rotation.  The  method  of  power  dissipation  could  be  the  result  of  either 
Coulomb  or  viscous  friction.  Viscous  friction  was  chosen  as  the  preferred 
method.  In  the  early  stage  of  the  design  it  was  thought  that  an  impeller 


Mt  = total  mass 

I = moment  of  inertia  about  axis  of  rotation 
z 

Ix  = moment  of  inertia  about  point  of  angular  deflection 
K&  = stiffness  refered  to  point  of  angular  deflection 
f = whirling  speed  ratio 
Vf8  * shaft  speed 
1 = shaft  length 

IP  = horsepower  required  to  produce  Torquewhlrl 
= damping  ratio 


could  be  placed  inside  of  a housing,  supported  on  the  same  shaft  that  was 
used  to  drive  the  impeller.  The  power  was  to  be  dissipated  by  pumping  the 
viscous  fluid.  There  were  numerous  problems  with  this  method,  particularly 
with  the  accurate  calculation  of  power  dissipation.  Consequently  this  idea 
was  dropped.  It  was  replaced  by  a system  based  on  thin  film  viscous  shear. 

In  this  thin  film  shear  idea,  the  housing  fit  very  close  to  a flat  rotating 
disk.  Power  would  be  dissipated  by  the  shearing  of  a very  thin  film  of  viscous 
fluid.  A good  deal  of  power  can  be  absorbed  by  this  technique.  At  the  same 
time  this  idea  was  being  investigated,  the  complex  interrelationships  of  the 
many  variable  design  parameters  were  not  understood,  and  the  method  appeared 
to  require  a system  of  impractical  size.  In  hindsight,  this  method  does  indeed 
seem  to  be  workable;  however,  a design  idea,  put  forward  at  a meeting*  in 
mid-January  of  '78,  forms  the  basis  of  the  current  system.  In  this  idea  the 
impeller  took  the  form  of  a conical  plug  (Fig.  A-3A).  The  plug  fit  inside  of 
a closely  mated  housing.  The  shape  of  the  outside  of  the  housing  was  not  known 
at  the  time,  but  was  assumed  to  be  boxlike.  Power  dissipation  in  this  system 
was  adjustable.  The  adjustment  could  be  effected  by  raising  or  lowering  the 
plug  thereby  changing  the  thickness  of  the  fluid  film.  Still,  at  the  time, 
the  size  of  the  actual  system  had  not  been  determined.  Mathematical  models  of 
different  sized  systems  were  developed  and  investigated.  The  relation  of  the 
various  design  variables  was  developed  from  these  models.  An  equation  relating 
several  of  the  parameters  was  used  to  generate  required  h.p.  dissipation  vs. 
speed  curves  for  the  various  models. 

For  all  calculations  the  fluid  viscosity  was  assumed  to  be  56  micro  reyns 
(Dow  Coming  silicone  1000).  The  film  thickness  was  taken  to  be  0.003  in.  and 


*The  meeting  was  comprised  of  Dr.  J.  M.  Vance,  Dr.  G.  N.  Sandor,  and  Mr.  Ken  Ard. 


CONICAL  PLUO  IN  MATING  HOUSING 
VISCOUS  SHEAB  DISSIPATION 


and  the  cone  plug  half  angle  was  30  degrees.  Throughout  the  entire  design 
process  the  whirling  speed  ratiw  was  limited  to  values  between  0.4  and  0.7, 
so  as  to  resemble  actual  field  measurements.  From  a review  of  the  resultant 
data,  one  system  was  chosen  as  the  most  desirable.  At  the  time  this  selection 
was  made  most  of  the  parameter  relations  were  not  known.  It  was  through  careful 
examination  of  the  effects  of  changing  parameters  that  the  relations  became 
apparent.  The  selected  system  was  continually  modified.  The  purpose  of  all 
modifications  was  to  increase  speed  and  to  reduce  required  horsepower.  The 
result  of  the  basic  modifications  was  an  overall  decrease  in  size  and  the  removal 
of  the  shaft  length  from  the  coupling  to  the  dissipation  housing.  A sketch 
of  the  basic  unit  as  perceived  at  this  time  is  shown  in  Fig.  A-4A. 

Problems  with  the  Basic  Unit 

For  any  machine,  the  percent  critical  damping  that  will  be  present  in  the 
final  construction  is  very  difficult  to  determine.  In  the  Torquewhirl  rig 
damping  is  very  important.  The  greater  the  damping,  the  greater  the  power 
required  to  produce  the  whirl.  Before  any  design  could  be  developed,  it  was 
necessary  to  find  some  approximation  of  the  percent  of  critical  damping  that 
would  be  present  in  the  basic  system.  To  gain  this  information  a simple  model 
was  constructed  and  test  measurements  were  made.  From  the  results  of  the  tests 
it  was  determined  that  the  damping  was  low  enough  to  result  in  an  acceptable 
power  requirement. 

The  system  in  its  present  form  will  dissipate  close  to  12  h.p.  The  imme- 
diate problem  is,  of  course,  how  to  get  rid  of  12  h.p.  worth  of  heat.  It  was 
first  hoped  that  the  thermal  capacity  of  the  housing  would  be  sufficient  to 
allow  the  time  necessary  to  take  data  without  excessive  temperature  rise.  This 
proved  not  to  be  the  case.  The  housing  would  reach  temperatures  in  excess  of 
250  degrees  F in  less  than  59  seconds.  When  the  effects  of  natural  convection 
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are  included,  no  significant  improvement  is  realized.  Forced  convection  was 
the  natural  alternative.  Fins  were  considered  necessary  in  order  to  achieve 
the  required  heat  transfer  rate.  Forced  flow  was  assumed  to  be  around  6 ft. /sec 
based  on  the  flow  velocity  of  a standard  blower.  Assuming  a steady  state 
temperature  of  197  degrees  F,  eight  fins  would  be  required.  The  final  design 
has  sixteen  fins  and  will  operate  well  below  175  degrees. 

Another  problem,  not  quite  as  obvious  as  the  heating  problem,  is  a result 
of  the  centrifugal  effect  imparted  to  the  fluid  by  the  rapidly  rotating  plug. 

This  centrifugal  effect  tends  to  drive  the  fluid  out  from  beneath  the  plug. 

The  result  of  this  motion  is  to  create  a suction  which  pulls  the  plug  down 
into  the  housing,  the  final  result  being  seizure.  This  problem  is  avoided  by 
drilling  return  ports  in  the  plug  itself  to  allow  the  pumped  fluid  to  return. 

The  pumped  fluid  would  be  slung  out  of  the  housing  if  a cover  plate  were  not 
included.  This  plate  serves  two  purposes;  first,  it  prevents  the  fluid  from 
being  thrown  out,  but  it  also  serves  as  a place  to  mount  the  upper  bearing  of 
the  cone-shaped  rotor.  The  upper  bearing  is  necessary  for  the  gap  control 
apparatus.  It  also  helps  ensure  proper  alignment  and  prevents  wobble. 

In  order  to  change  the  power  dissipation  independently  from  the  speed,  it 
is  necessary  to  change  the  fluid  film  thickness.  In  the  preliminary  design 
this  was  accomplished  through  the  tightening  of  a nut  below  the  housing  which 
raised  or  lowered  the  housing  with  respect  to  the  plug.  In  the  final  design 
this  same  basic  idea  is  applied  in  a different  form.  The  raising  or  lowering 
of  the  housing  is  done  by  using  two  wedges.  One  wedge  is  grounded  to  the  housing 
while  the  other  is  free  to  slide.  Due  to  the  angle  of  the  wedge  the  horizontal 
motion  is  converted  into  vertical  motion  of  the  housing.  A half  inch  of  wedge 
travel  is  required  to  produce  0.007  inches  of  film  gap  change. 


A coupling  selection  was  made  based  on  an  estimated  4 degrees  of  angular 
misalignment,  while  transmitting  70  ft.  lbs.  of  torque. 

Standard  universal  joints  could  not  be  used,  since  when  deflected  they 
do  not  transmit  at  constant  angular  velocity.  After  some  searching,  the 
Thomas  model  162  coupling  was  chosen.  It  is  of  the  standard  disk  pack  con- 
struction and  will  handle  the  expected  torque  and  misalignment. 

To  power  the  rig,  a 15  h.p.  aircraft  generator  was  purchased  to  serve  as 
a prime  mover.  The  basis  for  this  selection  was  primarily  price,  but  also 
because  some  experience  has  been  gained  with  this  type  of  power  drive. 

The  included  drawings  are  self-explanatory  and  provide  a detailed  de- 
scription of  all  the  parts  of  a test  rig  as  they  will  be  fabricated.  Some 
parts,  such  as  the  bearings  and  the  coupling,  have  either  been  ordered  or 
already  purchased.  The  major  portion  of  the  support  structure  and  safety  cage 
has  been  constructed. 

A general  idea  of  the  test  rig  can  be  derived  through  inspection  of  the 
assembly  drawings.  The  assembly  drawing  has  been  divided  into  two  parts  — 
the  upper  half  and  the  lower  half.  The  upper  half  assembly  shows  the  upper 
drive  shaft  (Dwg.  5),  the  upper  shaft  bearings,  and  mounting  plates  (Dwg.  3). 
The  mounting  plates  are  held  apart  by  spacer  blocks  (Dwg.  4) . The  upper  shaft 
is  rigidly  held  in  position  by  the  precision  Barden  bearings  and  the  all  steel 
support  structure  of  the  upper  half.  The  upper  mounting  plates  are  to  be 
positioned  during  construction.  They  will  then  be  firmly  clamped  to  each  other 
and  to  the  support  frame  by  four  one  half  inch  socket  headed  cap  screws  and 
nuts.  Further  positive  location  is  assured  through  the  use  of  four  one  quarter 
inch  dowel  pins  which  are  set  in  the  spacer  blocks  through  the  mounting  plates. 

The  upper  half  of  the  Thomas  coupling  is  attached  to  the  lower  end  of  the 
upper  drive  by  a key  and  setscrew.  A clear  representation  of  the  drive  shaft 


shaft  can  be  seen  in  Drawing  No.  5.  The  design  is  such  that  accurate  machining 
is  specified  only  where  necessary  along  the  shaft  (e.g.  bearing  seats).  The 
purpose  of  this  is  to  save  both  time  and  money. 

The  lower  portion  of  the  assembly  is  the  most  important  part  as  far  as 
Torquewhirl  considerations  are  concerned.  The  dynamic  parameters  of  this  half 
determine  the  power  and  speed  combinations  that  will  satisfy  the  Torquewhirl 
theory.  These  important  values  have  been  .carefully  calculated.  However,  the 
actual  parameters  are  expected  to  vary  slightly  from  the  calculations.  In  view 
of  this  fact,  some  items  necessary  for  rig  operation  have  not  been  included  in 
the  present  set  of  drawings.  The  actual  dimensions  of  these  small  and  inex- 
pensive parts  will  not  be  known  until  measurements  can  be  made  on  the  completed 
lower  assembly.  The  lower  assembly  (Dwg.  2)  consists  of  a finned  housing 
(Dwg.  10)  enclosing  a conical  plug  (Dwg.  9)  which  is  free  to  rotate  in  the 
housing.  There  is  a housing  cover  (Dwg.  8),  a slide  adjustment  assembly  (Dwg.  7), 
the  lower  drive  shaft  (Dwg.  12),  lower  shaft  bearings,  and  the  lower  half  of 
the  flexible  coupling  (Dwg.  6). 

The  lower  shaft  design  is  similar  to  that  of  the  upper  drive  shaft.  The 
reasons  for  this  type  of  design  are  again,  time  and  expense.  The  flexible 
coupling  is  pressed  onto  the  end  of  the  lower  shaft  and  keyed  to  resist  the 
load  torque.  Other  than  the  shaft,  bearings,  and  assorted  bolts,  the  entire 
lower  assembly  is  made  of  aluminum.  The  object  of  this  feature  is  to  lower 
the  weight,  and  therefore  moment  of  inertia,  of  the  lower  portion  of  the  system. 
This  lower  weight  is  required  in  order  to  keep  the  power  and  speed  combination 
required  to  produce  Torquewhirl  to  a reasonable  value.  (Note  Linking  in  Fig.  A-1A). 

The  conical  plug  is  also  keyed  to  the  shaft  in  order  to  stand  the  load 
torque.  The  plug  is  basically  a truncated  cone  of  30°  half  angle.  Holes  are 
drilled  into  and  through  the  cone,  to  lighten  it  and  to  provide  return  ports 
for  fluid  which  is  pumped  to  the  top  of  the  cone  due  to  its  rotation.  The  cone 
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fits  snugly  into  a mating  finned  housing  (Dwg.  10).  The  housing  holds  the 
shearing  fluid  and  is  finned  to  provide  for  adequate  cooling.  Cooling  of 
the  lower  housing  will  be  accomplished  by  the  use  of  a standard  blower  whose 
blast  will  be  directed  at  the  base  of  the  finned  housing.  There  is  a cap 
plate  (Dwg.  8)  that  covers  the  top  of  the  housing  to  prevent  spillage  of  the 
shearing  fluid.  This  cap  also  holds  the  upper  bearing  for  the  lower  shaft. 

There  is  an  assemblage  of  two  wedges  which  is  fastened  to  the  top  of  the 
housing  cap  plate  (Dwg.  7). 

These  wedges  are  used  to  provide  adjustment  of  the  film  thickness  of  the 
shearing  fluid.  This  adjustment  allows  the  independent  control  of  power  dis- 
sipation. Changes  in  the  film  thickness  are  effected  by  sliding  the  wedge  in 
a horizontal  direction.  The  wedge  angle  is  such  that  one  half  inch  in  wedge 
motion  results  in  seven  thousandths  of  an  inch  of  film  gap  change.  The  sliding 
wedge  is  held  in  position  by  an  overhanging  ledge  that  fits  over  the  fixed 
wedge.  The  adjusting  screws  also  aid  in  maintaining  the  wedge  position. 

A part  which  is  not  shown  in  this  assembly  drawing  is  the  drag  rod,  which 
is  used  in  adjusting  the  damping  value  of  the  whirling  motion.  The  drag  rod 
is  nothing  more  than  a solid  circular  rod  extending  downward  from  the  bottom 
of  the  finned  housing.  The  damping  will  be  adjusted  by  raising  or  lowering 
the  level  of  viscous  fluid  in  a container  in  which  the  lower  part  of  the  drag 
rod  is  immersed.  The  drag  rod  is  to  be  attached  to  the  bottom  of  the  housing 
by  adhesive. 

The  test  rig  as  described  is  expected  to  operate  at  speeds  between  one 
and  two  thousand  revolutions  per  minute.  Power  dissipation  is  calculated  to 
be  between  eight  and  twelve  horsepower.  Theory  predicts  that,  with  the  above 
limits  on  speed  and  power,  there  will  exist  numerous  combinations  that  will 


produce  Torquewhirl  in  the  system.  The  amplitude  and  stability  of  this  whirl 


ASS 


motion  are  to  be  measured  using  either  capacitance  or  magnetic  transducer* 
which  are  available.  It  is  expected  that  this  rig,  as  designed,  will  answer 
the  need  for  formal  verification  of  the  Torquewhirl  theory. 
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Caption  for  Figure  A-1A  Page  A25 


The  general  method  used  to  size  the  Torquewhirl  test  rig  is  outlined 

in  figure  A-1A,  Page  A25.  First  it  was  necessary  to  choose  a size  and 

shape  for  the  rig.  The  required  values  I and  I were  calculated. 

JC  z 

K.  was  known  between  limits,  and  was  used  along  with  W to  force  the  system 
y s 

within  the  bounds  defined  by  f.  If  KQ  and  W§  turned  out  to  be  reasonable, 
the  system  then  seemed  to  be  within  design  limits.  Next  it  was  necessary 
to  test  to  see  if  the  horsepower  required  to  produce  Torquewhirl  was  less 
than  the  twelve  horsepower  limit.  If  the  system  required  less  than  twleve 
h.p.,  it  was  examined  to  determine  if  it  would  actually  dissipate  the 
required  amount.  The  lower  the  h.p.  requirement  the  better  since  numerous 
combinations  of  speed  and  power  that  satisfy  the  Torquewhirl  theory  would 
then  exist  below  the  twelve  h.p.  limit.  The  selection  of  the  final  design 
was  an  itterative  process,  which  was  followed  until  all  requirements  of  the 
theory  were  satisfied. 
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Theory  and  Description 
of  Test  Rig 


Figure  B-l  is  a representation  showing  that  as  the  radial  tip  clearance 

of  the  blades  of  a compressor  varies,  there  is  a corresponding  variation  in 
efficiency.  Large  radial  clearances  cause  low  efficiency,  which  in  turn 
requires  increased  work  input  in  the  form  of  large  blade  forces  to  pro- 
duce the  required  static-pressure  rise.  The  output  required  of  an  axial - 
flow  compressor  is  such  that  the  axial  pressure  rise  is  the  same  for  every 
portion  of  the  circumference. 

The  local  clearance  in  the  X-direction  is  shown  in  Figure  B-2.  The 
following  is  taken  from  the  Journal  of  Engineering  for  Power  (Ref.  1): 


"Assuming  rotor  and  stator  to  be  perfect  circles, 
rotor  displacement  x will  result  in  a circumfer- 
ential variation  of  blade-tip  clearance.  The 
local  clearance  at  position  6 is 

6*6  ♦ x cos  0,  where  6 ■ mean  clearance, 

m m 

"Thermodynamic  stage  efficiency  is  a function  of 
tip  clearance,  the  variation  being  proportional 
to  ratio  of  eccentricity  to  blade  or  bucket  height. 


— The  local  efficiency  is 

" ■ V1  * r00595 

__  • 

where  nm  is  the  average  or  mean  efficiency  corre- 
_ spending  to  the  average  or  mean  clearance,  6^. 

H is  the  axial  height  of  the  blade.  For  turbines. 


8 has  a value  of  about  1.0  to  1.5.  For  compressors 


0 also  has  similar  values. 


Multistage  axial  compressors  with  some  stages 
close  to  stall  may  have  substantially  higher 
values.  This  results  from  the  effect  of  tip 
clearance  on  stall  pressure  ratio." 

"For  a turbine  with  constant  speed  and 
constant  entrance  and  exit  conditions,  the 
output  torque  is  directly  proportional  to  the 
local  efficiency.  The  local  torque  per  radian, 
Tq,  is  then  related  to  the  average  or  mean 
torqile,  T , by  the  relation 

Te  ■ V1  * r”*8)- 

The  component  of  force  taken  in  the  direction 
perpendicular  to  the  instantaneous  displacement, 
x,  is  per  radian 

AF  — HcosS 
y r 


The  net  force  on  the  rotor  perpendicular  to  the 
instantaneous  displacement  is 

Fy  ■ -f/2lT(l  ♦ ^os9)cos6de 


Total  torque 

w”  ve 


J‘\a  * jb-we 


(48) 
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FIGURE  B:l., HP..  VS,  UP  CLEARANCE 


FIGURE  8-.Z  r..  RESULT  AMT  FORCE  DUE  TO  BLADE 

ECCENTRICITY 


Combining  equations  (48)  and  (49)  gives 


F . Ifii.  I§_ 

y 2Hr  DpH 
where 


qx 


(50) 


The  total  torque  is  related  to  hp  and  engine 
rpm  by 


hp  x 33,000  x 12 
T ' f Zffrpm 


’ 63’000  {Si*  in‘lb 

Therefore 

. 63, 000 (hp) 8 
q Dp  x H x (rpm) 


Rotation  is  assumed  counterclockwise.  A 
positive  displacement  along  the  x-axis  causes 
a positive  force  in  the  y-direction.  A 
positive  displacement  along  the  y-axis  causes 
a negative  force  in  the  x-direction.  " 

Two  test  rigs  have  been  built  for  the  purpose  of  verifying  equation 
50.  The  first  test  rig  had  several  shortcomings  which  were  corrected  in 
the  second  test  rig  as  follows.  The  design  of  the  second  test  rig  incor- 
porated the  salient  features  of  the  first,  but  had  many  added  features, 
such  as  a heavy  rigid  frame  to  support  the  electric  motor,  shaft  and  rotor. 
The  new  frame  kept  the  vibration  down  to  a minimum.  A better  lock  mechanism 
was  provided  to  keep  the  shroud  surrounding  the  rotor  in  place  after  the 
eccentricity  was  introduced.  Sturdy  probe  mounts  were  added,  again  to 
keep  the  vibration  down  to  a minimum.  Slotted  motor  mounts  were  also  added 
to  give  the  motor  5 degrees  of  freedom  (3  translations  and  2 rotations  about 


horizontal  axes)  for  easy  aligning  of  the  motor  shaft,  quill  shaft  and 
rotor  shaft  with  each  other  and  for  the  use  of  different  length  quill 
shafts  (see  Sheet  No.  20).  Most  important,  what  prompted  the  evolution 
of  the  second  test  rig,  the  mean  clearance  between  the  blade  tips  and  the 
shroud  was  reduced  by  0.28  inches  to  a mean  clearance  of  O.OSO  inches.  It 
was  felt  that  the  large  clearance  caused  the  rotor  to  surge,  which  in 
turn  made  it  difficult  to  observe  any  aerodynamic  changes  due  to  Alford's 
force.  The  pitfalls  of  the  first  test  rig  were  thus  eliminated  in  the 
second  test  rig. 

The  eccentricity  was  imparted  to  the  rotor  by  moving  the  shroud  in 
a horizontal  direction.  It  varied  from  0.001  to  0.099  inches  and  could  be 
held  in  place  anywhere  in  between.  A probe  used  was  a Bently  proximeter 
probe  which  measured  horizontal  displacement  perpendicular  to  the  eccentricity. 
The  "balancing  top"  (see  Sheet  No.  16)  served  as  a surface  for  the  probe  to 
act. 

The  rotor  itself  was  mounted  on  an  elastic  beam  (see  Sheet  No.  18) 
which  allowed  it  to  move  in  a horizontal  direction  perpendicular  to  the 
blade  eccentricity  when  a small  force  was  applied  at  the  rotor  blades. 

Due  to  the  non-concentric  circumference  of  the  balancing  top,  it  was 
necessary  to  build  an  L.  C.  Low-Pass  filter  to  attenuate  the  A.  C.  component 
of  the  signal  produced  by  the  Bently  proximeter  probe.  The  signal  was  dis- 
played on  an  oscilloscope.  The  L.  C.  filter  consisted  of  a 10  henry  choke 
in  series  with  a lead  to  the  oscilloscope  and  a 50  microfarad  non-polarized 
capacitor  connected  across  the  2 leads  entering  the  oscilloscope  (in  a parallel 
configuration) . 


In  order  to  verify  equation  50  (Alford's  equation)  it  was  necessary  to 
know  the  hp  output  of  the  electric  motor  at  various  rpm-s.  A dynamometer  was 
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designed  and  built  for  this  purpose.  It  was  an  Eddy-Current  dynamometer  (Ref.  2) 
and  the  details  of  the  design  can  be  seen  on  Sheet  No.  22  of  the  assembly 
drawing.  The  eddy  current  brake  has  many  advantages  over  the  band  or  rope 
methods  of  motor-testing.  Among  these  are: 

a.  Sensitivity  combined  with  convenience  of  adjustment  (a  virtually 
infinite  adjustment  in  torque) 

b.  Exact  constancy  after  final  temperature  conditions  have  been 
reached 

c.  Heat  produced  in  energy-absorbing  element  not  conducted  to  the 
bearings 

d.  Absence  of  wear 

The  first  three  arise  from  the  nature  of  the  retarding  force  employed, 
and  from  the  simplicity  and  fineness  of  rheostatic  control.  The  continual 
stream  of  air  outward  past  the  copper  disc  carries  away  the  heat.  This  is 
the  advantage  referred  to  in  c. 

The  next  phase  of  the  research  on  this  project  will  be  to  verify  equation 
50.  The  previous  design  work  mentioned  above  will  lead  to  an  auspicious 
beginning  of  this  phase.  The  Alford's  equation  will  be  compared  with  the 
measured  plot  of  Alford's  force  vs.  motor  rpm  (eccentricity  held  constant) 
and  also  by  a plot  of  the  measured  Alford's  force  vs.  eccentricity  with  motor 
rpm  held  constant.  If  Alford's  equation  does  not  predict  the  actual  values, 
than  a new  equation  will  be  derived  using  the  experimental  data  obtained  from 
the  Alford  test  rig. 
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Appendix  C 

Torquewhirl  of  a Flexible  Shaft  Rotor 

Reference  1 gives  an  exact  functional  solution  to  the  nonlinear  dif- 
ferential equations  of  motion  for  a torque -loaded  disk  mounted  on  a straight 
rigid  shaft,  overhung  from  a flexible  joint.  This  solution  shows  that  load 
torque  alone  can  produce  nonsynchronous  whirling,  with  amplitudes  dependent 
on  the  ratio  of  torque  to  damping. 

A model  which  more  closely  approximates  a real  machine,  with  a single 
overhung  disk,  is  shown  in  Figure  C-l.  The  shaft-  of  length  1 is  assumed 
to  be  flexible,  but  its  mass  is  neglected  as  being  insignificant  conpared 
to  the  disk. 

Whereas  the  model  of  reference  1 has  three  degrees  of  freedom,  the 
model  of  Figure  C-l  has  five.  If  a driving  torque  along  the  Z axis  is 
reacted  by  a load  torque  aligned  with  the  disk  axis,  then  it  might  be 
supposed  that  the  flexible  shaft  model  has  a solution  for  nonsynchronous 
whirl  similar  to  that  of  reference  1. 

To  investigate  this,  consider  generalized  coordinates,  x,  y,  e,  ♦, 
(Figures  C-l  and  C-2).  The  advantage  of  the  Euler  angles  is  that  whirling 
is  represented  by  constants , simplifying  the  search  for  a solution.  The 
x,y  coordinates  do  not  have  this  advantage,  and  a whirling  solution  requires 
harmonic  (time-dependent)  functions.  Therefore,  consider  replacing  x,y  by 
R,a,  where  R is  the  projection  of  CS  on  the  xy  plane  and  a is  the  angle 
XOG.  the  conplete  set  of  coordinates  is  then  R,  a,  e,  ♦,  i|»,  and  the  kinetic 
energy  is  expressed  as 


where 


The  potential  energy  is  elastic  (strain)  energy  stored  in  the  shaft 
and  will  be  a function 


The  potential  energy  is  more  easily  written  as  V(X,  Y,  3 , a , «j/) 

x y 

where  a , a are  disk  rotations  about  x and  y.  Therefore,  it  may  be 

A / 

productive  to  derive  transformations  from  x,  y,  a , a to  R,  a,  0, 

A 7 

Since  <|;  is  pure  twist,  it  may  appear  in  an  uncoupled  term. 

The  virtual  work  of  the  load  torque  produces  the  same  generalized  torques 
and  as  in  reference  1.  Most  of  the  external  damping,  however,  will 
act  on  ft  and  Rd,  rather  than  on  6 and 


Considerable  simplification  is  acccnplished  by  assuming  that  the  shaft 
remains  in  the  plane  BOG  when  it  bends,  thus  reducing  the  degrees  of  freedom 
to  four.  An  attractive  set  of  coordinates  is  0,  $,  * (see  Figures  C-l  and 
C-2)  also,  let  * ■ $ ♦ ir/2  and  the  following  relationships  can  be  written: 
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tan«  - - , 


(C-3) 


x^  ♦ y 


(C-4) 


w ■ ♦ ♦ ♦ 


(C-5) 


The  body  rates  are: 


wx  * 4 sinesiity  + 6 cos<|>, 


(C-6) 


Wy.  * i sinecosij)  - 6 simp, 


IC-7) 


w ■ $ + $ cose. 

Z 


(C-8) 


Note  that  w f w2< 


The  kinetic  energy  can  now  be  written  as 


T - sine  + -tg  6 cose)2  + ($  te  sine)2]  + ^ Ix[e2  + $>2  sin2e] 


\ I2(^  + 4,  cose)' 


(C-9) 


for  a symnetric  rotor  (Ix  * 1^). 


The  potential  energy  is 


V - V(e,  lQ) 


(C-10) 


The  virtual  work  of  the  damping  forces  is 


6w  > -c.tt  sine  + cose](6£  sine  + l 6 e cose) 


' Cd^e  sin0](£e  sin0)  5*» 


and  the  virtual  work  of  the  driving  torque  and  load  torque  is  the  same  as  in 


reference  1. 


... 


The  three  differential  equations  of  motion  are  obtained  by  substituting 


L ■ T - V and  the  Qi  into 

where  qx  -*e,  q2  - e,  q3  - ♦ , q4  • *, 

The  solution  should  be  similar  to  that  of  reference  1. 
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Appendix  D 

Torquewhirl  Stability  Predictions  for  Multi-Disk  Rotors 
by  the  Transfer  Matrix  Method 

Lund  (reference  14)  has  extended  the  Myklestad-Prohl  transfer  matrix 
method  to  rotor-bearing  systems  which  include  damping  and  destabilizing  cross - 
coupled  stiffness  and  danping  coefficients.  Whereas  the  Myklestad-Prohl  method 
yields  only  the  imaginary  part  of  the  eigenvalues  (i.e.  the  natural  frequencies), 
the  Lund  method  yields  complex  eigenvalues  (i.e.  both  the  natural  frequencies 
and  the  logarithmic  decrement,  which  is  a stability  predictor). 

Currently  the  greatest  limitation  of  Lund's  method  (as  with  all  other 
stability  analyses)  is  the  lack  of  accurate  information  about  the  types  of 
destabilizing  excitations  which  exist  in  real  machines  and  which  therefore  are 
to  be  used  as  input  to  the  computer  program.  Torquewhirl  (reference  1)  has 
recently  been  identified  as  one  of  these  excitations.  To  put  the  torquewhirl 
forces  into  a Lund  stability  analysis,  they  must  be  formulated  in  terms  of 
stiffness  and  damping  coefficients. 

The  coordinates  used  in  Lund's  analysis  are  shown  in  Figure  D-l.  In  gen- 
eral, the  stiffness  and  danping  coefficients  are  the  matrices  which  define  the 
forces  and  moments  on  each  disk  in  the  x and  y directions. 

For  example,  the  force  on  a disk  in  the  x direction,  due  to  disk  displace- 
ment and  velocity,  is 

Fx  - - v • V - V - V • Kx«»  - v - ^ • V 


and  the  matrix  equation  for  the  forces  in  all  directions  (on  a single  disk) 

is 
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where  Fq  and  Fg  are  actually  moments  on  the  disk. 

It  can  be  seen  that  there  are  a total  of  thirty  two  stiffness  and  damping 
coefficients  defining  the  forces  and  moments  on  each  disk.  The  off-diagonal 
elements  are  called  the  cross -coupled  coefficients. 

Fot  the  rigid-shaft,  flexible  joint,  torquewhirl  model  in  reference  1, 
the  disk  rotation  and  translation  coordinates  are  not  independent  (i.e.  they 

are  related  by  a kinematic  constraint).  The  constraint  equations  are  R - te 

2 2 2 

in  polar  coordinates  or  x ♦ y ■ (l  sine)  in  inertial  coordinates. 

Since  the  equations  and  generalized  forces  in  reference  1 are  written  in 
terms  of  Euler  angles  + , e,  and  ♦ , a coordinate  transformation  is  required 
to  derive  the  stiffness  and  danping  coefficients,  as  follows  (see  Figure  P-2): 


First,  the  generalized  torque  Q.  of  reference  1 can  be  expressed  as  a 

♦ 

tangential  force  (see  Figure  D-2) 


% - ty*. 

and  the  generalized  moment  Q.  can  be  expressed  as  a radial  force 
Fr  ■ Qe/£  cose. 

It  is  (or  F^)  which  contains  the  destabilizing  torquewhirl  forces. 
The  transformation  to  x and  y is 


04 


I X 

I I 

J.. 

1 

I I 

I J. 


where 


Fx  ■ Fr  cos*  - F^  sin* 


F - Fr  sin*  ♦ F^  cos* 


cos*  - 


(x‘  ♦ y ) 


(x"  ♦ y ) 


For  the  aerodynamic  case  (see  reference  1) , the  generalized  torque  Q 


Q*  - Ts  - C^(*  + * cose)"2  cose  - Cd(£3  sin3e)  *2, 


where 


Ts  - shaft  torque 

^ ■ disk  load  coefficient 

Cj  - nonlinear  damping  coefficient 

The  shaft  speed  w is 


wg  - * ♦ *. 


The  shaft  torque  equals  the  disk  load  torque  so  that 

*) 

Tg  • (*  + * cose) 

Therefore  the  generalized  torque  Q . can  be  expressed  as 

* 

- T$  [1  - cose]  - Cd  {f}  sin3e)  *2  . 


Rigid-Shaft  Model  (Reference  1) 

The  destabilizing  part  of  is  the  first  term  from  the  expression  just 

above.  The  equivalent  tangential  force  is 
Q.  T (1  - cose) 


The  required  relationships  between  R,  e,  and  a,  e,  and  x,  y are 


Therefore,  in  terms  of  x and  y,  we  can  write 


Keeping  only  the  first  two  terms  of  the  cosine  series  yields 


By  inspection,  it  can  be  seen  that  the  cross-coupled  stiffnesses  are 


It  is  interesting  to  note  that  and  have  the  sane  form  as  Alford’s 
coefficients  for  the  effect  of  tip  clearance  asymnetry  in  axial  flow  turbo- 
machinery (reference  7). 

As  a crude  approximation,  the  above  coefficients  (K^,  K ) could  be  used 
for  a flexible- shaft  model  by  taking  l as  the  axial  distance  from  the  disk 
plane  to  a virtual  pivot  point,  as  determined  by  the  mode  shape. 


Flexible- Shaft  Model 


When  the  rotor  flexibility  is  distributed  along  the  shaft,  there  is  no 
constraint  between  x,  y and  a,  6 (see  Figure  D-l) . Therefore  the  virtual 
work  of  the  destabilizing  torque  must  be  written  in  terms  of  x,  y,  a,  and  6 
This  requires  that  the  virtual  displacement  6$  be  written  in  terms  of  6x, 


Since  the  rotations  of  a disk  are  independent  of  the  translations,  it 
can  be  seen  that  6x  and  6y  will  not  appear  in  6$.  This  leaves  only  the  vir 
tiial  angular  displacements  5a  and  6$  to  be  considered.  A kinematical  rela- 
tionship between  ♦,  a,  and  8 is  thus  required. 

Assuming  this  can  be  obtained,  we  proceed  as  follows: 


For  a multi-disk  rotor,  the  destabilizing  torque  on  the  n disk  is 


where 


fraction  of  the  shaft  torque  accounted  for  by  the  nth  disk 
(a^  + (see  Figure  D-l) 


Therefore  the  virtual  work  of  the  shaft  torque  on  the  n disk  is 


r. 

i. 

j. 

l 
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Derivation  of  the  Characteristic  Polynomial  from  the  Transfer  Matrices 

In  the  computer  analysis  of  reference  14,  the  complex  roots  of  the 
characteristic  equation  are  obtained  by  a numerical  iteration  algorithm. 

In  some  cases  involving  a large  number  of  mass  stations  along  the  rotor, 
the  algorithm  has  been  known  to  miss  certain  roots.  Also,  the  roots  are  not 
determined  in  any  particular  order,  and  the  higher  order  roots  are  less 
accurately  determined. 

These  shortcomings  suggest  that  it  would  be  advantageous  to  have  the 
characteristic  polynomial  available,  so  that  more  efficient  computer  sub- 
routines ("rootfinders")  could  be  used  to  determine  the  roots.  Of  course, 
the  polynomial  could  be  obtained  by  expanding  the  determinant  of  the  large 
eigenmatrix,  but  this  would  require  a very  large  computer  storage  array, 
thus  losing  one  of  the  advantages  of  the  transfer  matrix  method. 

Also,  electrical  engineering  and  controls  technology  have  produced  a 
wide  and  powerful  array  of  stability  investigation  techniques  which  use  the 
characteristic  polynomial  as  a starting  point.  There  is  no  reason  to  doubt 
that  many  of  these  techniques  would  be  useful  for  rotor  dynamic  stability 
studies. 

Derivation  of  the  characteristic  polynomial  by  multiplication  of  fifteen 
or  twenty  (four  by  foui)  transfer  matrices  for  a rotor  is  a challenging  task. 
Attempts  to  derive  functional  expressions  by  direct  multiplication  have  been 
successful  only  for  four  or  less  mass  stations  an  the  rotor.  A numerical 
scheme  is  outlined  in  reference  15,  but  is  implemented  only  for  a simple 
undamped  beam  vibrating  in  one  plane.  Therefore  the  most  recent  efforts  on 
this  Grant  have  been  directed  toward  developing  a formal  scheme  to  generate 
the  nunerical  values  of  the  polynomial  coefficients  on  a digital  computer. 
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To  develop  such  a scheme,  it  was  found  helpful  to  begin  with  a somewhat 
simpler,  although  completely  analogous,  problem.  This  is  the  generation  of 
polynomial  coefficients  for  torsional  vibration  of  a multi-disk  rotor. 

Figure  D-3  shows  the  model  used. 

The  transfer  equation  across  the  n^1  inertia  station  is 

T;  - Tn  * sVn  * SV„  * Vn 


where  S is  the  complex  eigenvalue. 

Across  the  n^  (massless)  shaft  segment,  the  equations  are 


TA 


and 


t; 


>'  . , a 


A numerical  confutation  of  the  eigenvalues,  by  HOlzer's  method,  starts 
with  an  assumption  for  S along  with  setting  the  arbitrary  amplitude  6^  ■ 1 
at  the  left  hand  end.  The  transfer  equations  above  are  then  used  repetitively, 
from  left  to  right  along  the  length  of  the  rotor,  to  calculate  the  right  end 
boundary  torque.  If  this  torque  is  correct,  the  assunption  for  S is  known  to 
be  an  eigenvalue. 

Here,  the  transfer  equations  are  to  be  used  to  derive  the  coefficients 
of  the  characteristic  polynomial,  which  has  the  form 


aQ  ♦ ajS  ♦ a2S^  ♦ a^S^  ♦ . . . ♦ ®n-2^  ^ 
for  an  N- inertia  rotor  with  zero  boundary  torques.  This  can  be  achieved  by 
either  of  the  following  two  methods.  (The  case  with  no  damping  is  illus- 
trated for  simplicity). 
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Method  1: 

The  amplitudes  and  are  eliminated  by  substitution  from  the  trans- 
fer equations  to  yield  (for  S - iw) 


TA 


Start  at  the  left  end  where  n ■ 1 and  work  to  the  right  by  substitution. 
For  exanqsle,  substitute  T'  for  T^2  and  T’  for  to  generate  Tj.  T|  and 
T2  are  made  especially  siaple  by  the  asstmption  0^-1.  They  are 


and 


In  the  computer  program,  all  parameters  except  w have  numerical  values. 
Coefficients  of  like  powers  of  w are  sunned  into  an  array  as  the  multiplica- 
tion progresses  to  the  right  end,  where  the  boundary  torque  is  zero. 


Method  2; 

The  polynomial  coefficients  of  the  four- inertia  rotor  can  be  obtained 
in  the  following  form,  which  suggests  a pattern  for  the  more  general  N- inertia 
system. 
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Extending  the  pattern  to  a five-inertia  rotor  yields 

1 j 

_ kik2kSk4^Il  + l2  * l3  + *4  + V 

“o 

WsVs 

jj 

etc. 

It  appears  that  methods  analogous  to  the  two  methods 

above  can  be 

Jj 

developed  to  generate  the  characteristic  polynomial  for  rotor  whirl  stability 

f 

predictions. 
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Plgure  D-lt  Virtual  Work  in  Terms  of  Coordinates 
Used  In  Lund’s  Program 


The  position  of  the  disk  Is  described 
by  x,y,c (z  Is  negligible  andu>is 
constant) 

x & y give  position  of  CO 

<*  kfi  give  angular  orientation  of  the 
disk  (Lund  oalls  themed^) 


